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Abstract 

Let each point of a homogeneous Poisson process on M indepen- 
dently be equipped with a random number of stubs (half-edges) ac- 
cording to a given probability distribution on the positive integers. 
We consider schemes based on Gale-Shapley stable marriage for per- 
fectly matching the stubs to obtain a simple graph with degree distri- 
bution 11. We prove results on the existence of an infinite component 
and on the length of the edges, with focus on the case ^({2}) = 1. In 
this case, for the random direction stable matching scheme introduced 
by Deijfen and Meester we prove that there is no infinite component, 
while for the stable matching of Deijfen, Haggstrom and Holroyd we 
prove that existence of an infinite component follows from a certain 
statement involving a finite interval, which is overwhelmingly sup- 
ported by simulation evidence. 



1 Introduction 

Let V he a homogeneous Poisson process with intensity 1 on R'^ and fi a prob- 
ability measure on the strictly positive integers. We shall study translation- 
invariant simple random graphs whose vertices are the points of V and where, 
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Figure 1: The stable multi-matching (top) and the random direction stable 
multi-matching, for 10 vertices on a finite interval, with 2 stubs per vertex. 

conditional on P, the degrees of the vertices are i.i.d. with law /x. Previously, 
Deijfen [2] has studied achievable moment properties for the edges, and Dei- 
jfen, Haggstrom and Holroyd [5] have studied the question of whether the 
graph contains a component with infinitely many vertices. In the latter work 
a particular matching scheme, called the stable multi-matching, was intro- 
duced, leading to a number of challenging open questions. Here we restrict 
to d = 1 and the focus is on the case /i({2}) = 1, one of the simplest cases for 
which the question of existence of an infinite component is non-trivial. For 
the stable multi-matching and a variant of it with prescribed random stub 
directions, we prove results on the component structure and on the length 
of the edges. Figure 1 shows schematic pictures of the two matchings, which 
are described below. 

First we formally describe the objects that we will work with. Write 
[V] := {x G M : V{{x}) > 0} for the support, or point-set, of V. Let 
be a random integer-valued measure on M with the same support as V, and 
which, conditional on V, assigns i.i.d. values with law jj, to the elements of 
[P]. The pair (V,C,) is a marked point process with positive integer- valued 
marks. For x G [V] we write for C,{{x}) and interpret this as the number 
of stubs at vertex x. 

A matching scheme for a marked process (V,^) is a point process Ai 
on the space of unordered pairs of points in M, with the property that al- 
most surely for every pair {x,y) G [Ai] we have x,y E [P], and such that 
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in the graph G = G{V,M.) with vertex set \P] and edge set [M], each ver- 
tex X has degree D^. Our primary interest is in the connected components 
of G. The matching schemes under consideration will always be simple, 
meaning that G has almost surely no self-loops and no multiple edges, and 
translation-invariant, meaning that M. is invariant in law under the action 
of all translations of M. We say that a translation-invariant matching is a 
factor if it is a deterministic function of the Poisson process V and the mark 
process ^, that is, if it does not involve any additional randomness. We write 
P and E for probability and expectation on the probability space supporting 
the random triplet ("P,^, A^). 

Let (P*, M*) be the Palm versions of (P, ^, M) with respect to V and 
write P* and E* for the associated probability law and expectation operator. 
Informally speaking, P* describes the conditional law of ("P, M) given that 
there is a point at the origin, with the mark process and the matching scheme 
taken as stationary background; see e.g. [HI Chapter 11] for more details. 

Since P is a Possion process, we have [V*] = [V] U {0}. 

We now define the two matching schemes that will be analyzed in the 
paper. 

Stable mult i- matching 

The concept of stable matching was introduced by Gale and Shapley [5]. 
It has been studied in [7] and [8] in the context of spatial point processes 
(with //({I}) = 1 in our notation). A natural generalization to other degree 
distributions /i was introduced in [3] and is referred to as the stable multi- 
matching. Formally, a matching scheme M. is said to be a stable multi- 
matching if a.s., for any two distinct points x,y E [V], either they are linked 
by an edge or at least one of them has no incident edges longer than \x — y\. 
Here and throughout, distance and edge length refer to the Euclidean norm 
I ■ I on M. 

We will restrict our attention to the case when P is a Poisson process. 
For this case, it was proved in [3l Proposition 2.2] that there is an a.s. unique 
stable multi-matching, which moreover can be constructed by the following 
iterative procedure. First connect all mutually closest pairs of points in 
[P] and remove one stub from each of these point. Then call two points 
compatible if they do not already have an edge between them and if both 
of them have at least one stub left. Connect all mutually closest compatible 
pairs and remove one stub from each of the points just matched. Repeat 
indefinitely. See [H Propostion 2.2]. 
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Random direction stable multi-matching 

We introduce a variant of stable multi-matching where the directions of the 
edges are prescribed independently of the Poisson process. As described 
above, the process ^ assigns a mark Dx to each point x G \P]. Let if) be 
a second mark process which, conditionally on V and assigns an integer 
Rx ~ Binomial (Z^a;, 1/2) independently to each point x E \P]. We think of 
Rx as the number of stubs incident with x that are to be matched to the right 
of X. If X < y, and is an edge of a matching scheme M., we call (x,?/) 

a right-edge of x, and a left-edge of y. A matching scheme Ai is now said 
to be a random direction stable multi-matching if each point x G [P] 
has exactly Rx incident right-edges and if a.s., for any two distinct points 
x,y E [V] with X < y, either they are linked by an edge, or x has no incident 
right-edges longer than \x — y\, ot y has no incident left-edges longer than 
\x - y\. 

Let each point x G [V] be equipped with Rx stubs pointed to the right and 
Lx := Dx — Rx stubs pointed to the left, and consider the following iterative 
procedure for matching right-stubs to left-stubs. First consider all pairs of 
consecutive points in [P]. Create an edge between every such pair x < y such 
that X has at least one right-stub and y has at least one left-stub, and remove 
the corresponding stubs. Then consider pairs of points in [V] with precisely 
one point in [P] in between them. Create an edge between every such pair of 
points X < y such that x has at least one right-stub and y has at least one left- 
stub left, and remove the corresponding stubs. Continue indefinitely, with 
pairs of points separated by an increasing number of points. This procedure 
has previously been studied in We show in Section [2] that it leads to the 
unique stable multi-matching subject to the prescribed (random) directions 
for the edges. 

1.1 Results 

In this section we collect the main results. The proofs are then given in 
Section [31 The first result concerns uniqueness of the infinite component. 

Proposition 1.1. For a Poisson process on M and any degree distribution, 
in the stable multi-matching and the random direction stable multi-matching, 
there is at most one infinite component. 

The next result asserts that, in the case /i({2}) = 1 of two stubs per 
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vertex, the random direction stable multi-matching has no infinite compo- 
nents. For other degree distributions the existence of an infinite component 
remains an open question. Part (b) of the theorem however provides some 
information on the edge length. See [H Theorem 4.1] and [71 Theorem 2] for 
related results. 

For X G [V], let denote the average length of all edges incident to 
X, and write X = Xq for the value at the origin in the Palm version of the 
process. 

Theorem 1.2. For a Poisson process on M, consider the random direction 
stable multi-matching. 

(i) For /i({2}) = 1, almost surely there is no infinite component. 

(ii) For any degree distribution with bounded support, we have 
E*[Xi/2] = oo. 

Turning to the stable multi- matching, it was proved in [3l Theorem 1.2(b)] 
that there is no infinite component when the only possible values for the 
degrees are 1 and 2, with a strictly positive probability of degree 1. In 
d > 2 ii was also proved that there is an integer k = k{d) such that if all 
vertices almost surely have degree at least k, then there is almost surely 
an infinite component, [3l Theorem 1.2(a)]. Note that, by ergodicity, the 
event that there exists an infinite component has probability or 1 for any 
degree distribution. The following result relates the existence of an infinite 
component for the case /i({l, 2}) = 1 in = 1 to a certain property concerned 
with the lengths of the edges. Let denote the length of the longest edge 
incident to x E [P], say that x desires a site y G M if |?/ — a;| < and write 
for the number of points in [P] that desire the origin. 

Theorem 1.3. For a Poisson process on M, consider the stable multi- 
matching. 

(i) For any degree distribution, if there is no infinite component, then 
N = oo almost surely. 

(ii) ///i({l,2}) = 1 and there is an infinite component, then N < oo almost 
surely. 

For degree distributions with /i({l,2}) = 1, existence of an infinite compo- 
nent in the stable multi- matching is hence equivalent to < oo. On the 
other hand, is related to edge lengths, as follows. Write M = Mq. 
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Lemma 1.4. For any translation-invariant matching scheme, we have that 
E*[M] <oo if and only ifE[N] < oo. 

In view of this relation, £'*[M] < oo would imply that < oo, and 
thereby establish the existence of an infinite component for /_i({2}) = 1 in 
the stable multi-matching. However, the best result we have in this direction 
is the following, which applies in any dimension d > 1 (the stable multi- 
matching is defined analogously in all dimensions; see [3]). 

Proposition 1.5. For a Poisson process of intensity 1 on W^, and any de- 
gree distribution with hounded support, in the stable multi-matching we have 
P*(M > t) < ct~'^ for some c G (0, oo) (depending only on d and the bound 
on degree). 

A "statistical proof" of percolation. It is not rigorously known whether 
the stable mult i- matching with /i({2}) = 1 will have an infinite component 
in d = 1. However, in Section S] we present compelling evidence that this is 
indeed the case. Specifically, we will define a certain event Gl in terms of a 
Poisson process on the bounded interval [0, L]. We will prove rigorously that 
for any L > 0, 

¥{Gl) > 0.968 implies existence of an infinite component. 

On the other hand, since Gl is defined in terms of a bounded interval, its 
probability can be estimated by Monte-Carlo simulation. Such simulations 
provide overwhelming statistical evidence that 

PlGisooo) > 0.968. 

The random direction stable multi-matching and the stable multi-matching 
are hence qualitatively different: when the directions of the stubs are pre- 
scribed randomly, there is no infinite component, while when the directions 
are prescribed by the positions of the Poisson points (as in the stable multi- 
matching) there is an infinite component. Figure 2 shows simulation pictures 
of the random direction stable multi- matching and the stable multi- matching, 
respectively, with fi{{2}) = 1 in = 1. 

Say that an edge (x, y) G [M.] crosses a site z^WifxKzKy. Our 
last result is the following. 
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Figure 2: The stable multi-matching (top) and the random direction stable 
multi-matching, with 2 stubs per vertex, for 500 uniformly random points on 
an interval. 

Proposition 1.6. For a Poisson process on M, and any degree distribution 
fi whose support includes some odd integer, for any factor matching scheme, 
the number of edges that cross the origin is infinite. 

If the number of edges that cross the origin if infinite, then clearly also 
= oo. Hence, appealing to Lemma II. 4[ Proposition 11.61 implies that 
E*[M] = oo in any factor matching scheme for degree distributions whose 
support contains an odd integer. Since the stable multi-matching is a factor, 
combining Proposition II. 61 with Theorem ll.3( b) gives an alternative proof (in 
d = 1) oi the result of [3] that the stable multi-matching does not percolate 
when the only possible values for the degrees are 1 and 2 and the probabil- 
ity of degree 1 is strictly positive. For degree distributions with support on 
larger values this approach is inconclusive, since Theorem 11.3( b) does not 
apply. 

The rest of the paper is organized as follows. In Section 121 a few prelim- 
inary results are collected. The above results are then proved in Section |3l 
In Section HI existence of a infinite component in the degree 2 case in d = 1 
is shown to follow from the assertion that a certain finite event has large 



enough probability, which is convincingly supported by Monte- Carlo simu- 
lation. Section |S] contains examples demonstrating that there is no general 
relation between the edge length and the existence of an infinite component 
valid for any matching scheme, and that for point processes other than the 
Poisson process, both percolation and non-percolation are possible for the 
stable multi-matching in the degree 2 case. Finally, in Section E] some di- 
rections for further work are presented. For background on the problem we 
refer to [3], Section 2.1] and [21 Section 1]. 

2 Preliminaries 

We first show that the iterative procedure described for the stable multi- 
matching with random directions leads to the unique stable multi-matching 
with the prescribed directions for the edges. 

Proposition 2.1. Let (V,^,ip) be a doubly marked Poisson process. Almost 
surely, the iterative multi-matching procedure described in the introduction 
exhausts the set of stubs, and the limiting graph (after an infinite number 
of iterations) is a random direction stable multi-matching. No other such 
matching scheme exists. 

Proof. Let P', (respectively P/) be the process of points with at least one 
unmatched right-stub (left-stub) on them after the above matching procedure 
is completed. By symmetry and V'l have the same intensity and they are 
both ergodic point processes. Hence either both have a.s. infinitely many 
points or both have a.s. no points. The first option however would produce 
a contradiction, since the iterative procedure could then be applied to the 
remaining configuration of stubs giving rise to edges that would have been 
created already in the original procedure. 

That the resulting multi-matching is stable subject to the prescribed (ran- 
dom) directions follows from the definition: an unstable pair of points - that 
is, a pair x and y with x < y with no edge between them and where x 
(y) has an edge connected to the right (left) of y (x) - would have had an 
edge created between them at some stage of the matching procedure. That 
it is the unique matching with this property follows by induction over the 
stages in the algorithm to show that each edge that is present in the resulting 
configuration must be present in any stable matching of the stubs with the 
prescribed directions. □ 
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We proceed by formulating a version of the mass transport principle suit- 
able for our needs. For background, see [T]. A mass transport is a random 
measure T on (M"^)^ that is invariant in law under translations of M.'^. For 
Borel sets A,B G M*^, we interpret T{A, B) as the amount of mass trans- 
ported from A to B. Write Q for the unit cube [0, 1)'^. 

Lemma 2.2 (Mass Transport Principle). Let T be a mass transport. Then 

ET{Q,R'^) =ET(M'^,g). 

Proof. 

ET{Q, R'^) = J2 ET(Q, Q + z) = ^^{Q - z,Q) = ET(M'^, Q) . □ 

Lemma [1.41 is now easily established using the mass transport principle. 

Proof of Lemma \1.4\ Consider the mass transport in which each point x G 
[V] sends out mass 2M^, and distributes it uniformly to the interval (x — 
Mx,x + Mj,). The expected mass sent out from the unit interval [0,1) 
equals 2E*[M]. On the other hand, the mass received by the unit inter- 
val is Nx dx, where A^^^ denotes the number of points that desire x G M. 
Hence the expected mass received by the unit interval is EA^. The result 
hence follows from the mass transport principle. □ 

Next we observe that an infinite component in a translation-invariant 
matching scheme must be unbounded both to the right and to the left, that 
is, for any r G M"^ it must contain points both to the right of r and to the 
left of — r. 

Lemma 2.3. A translation-invariant matching scheme almost surely cannot 
have an infinite component that is unbounded in only one direction. 

Proof. Assume that there is a matching scheme that with positive probability 
gives rise to an infinite component that is unbounded in only one direction, 
say to the left, and consider the mass transport in which each vertex in such 
an infinite component sends mass 1 to the rightmost point in the component. 
With positive probability such a rightmost point is located in the unit inter- 
val, which then receives infinite mass. The expected mass sent out from the 
unit interval is however bounded by 1, so we have a contradiction with the 
mass transport principle. □ 
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Finally, the following result will be of use in proving Theorem 11.3( b). 



Lemma 2.4. Let T be a translation-invariant simple point process of finite 
intensity on M. For x G [T], write for the maximum of the distances from 
X to the nearest point of [T] on the left and the nearest point on the right. 
The number of points x G [F] with > \x\ is finite almost surely. 

Proof. Consider the mass transport in which an interval sends out mass equal 
to its length, and the mass sent out by an interval (x, y) between consecutive 
points X < y of [F] is distributed uniformly to the interval {x — {y—x),y + {y — 
x)). If there were infinitely many points x G [F] with Zx > \x\, then the unit 
interval would receive infinite mass, which conflicts with the mass transport 
principle, since the mass sent out from the unit interval equals 1. □ 

3 Proofs 

We now proceed to prove the results in Section 11.11 starting with the unique- 
ness of an infinite component. We say that two edges (a, b) and (c, d) in [Ai] 
cross each other ii a < c < b < d. 

Proof of Proposition Observe that both matching schemes have the prop- 
erty that, 

if two edges (a, b) and (c, d) cross each other, then 
the edge (c, b) must also be present in the matching. 

- this follows from the definitions of the matching schemes. Two distinct 
components hence cannot have crossing edges. However, by Lemma 12. 3[ 
any infinite component must be unbounded in both directions. Hence two 
distinct infinite components would necessarily have crossing edges. □ 

Proof of Theorem \1.2\ (i). Recall that (respectively R^) is the number of 
edges incident with x G [P] that are connected to the left (right) of x. Let 
/i({2}) = 1. We call x a bird if L^. = i?^. = 1, a left-beak if L^, = 2 and a 
right-beak ii = 2 (see Figure 1). Let (■ ■ ■ <)xi <■■■< (<■■■) G [P] 
be the ordered vertices of some (finite or infinite) component of the stable 
multi-matching (recall by Lemma 12.31 that a component is either finite or 
unbounded in both directions). Clearly, if the component is finite, then its 
leftmost point is a right-beak and its rightmost point is a left-beak. We 
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claim that if ght-beak, and not one of these extreme points of the 

component, then Xj+i is a left-beak. To check this, let the two edges from 
Xi have their other endpoints at xj and Xk, where Xi < Xj < Xk- We claim 
that the other neighbour of xj must lie left of Xi. To see this, follow the 
path formed by the cluster starting with the edge {xi,Xj) - eventually we 
must leave the interval since the cluster contains points to the left 

of Xi- When we do so, it is via an edge that crosses {xi,Xk). Unless it is 
the first edge encountered after xj, this entails a violation of ([1]). Thus Xj 
is a left beak. Now there cannot be any further vertices of the cluster in 
the interval since such a vertex would have an incident edge crossing 

(^X{^Xj^^ again contradicting ([T]). 

Therefore, the non-extreme vertices of a component consist of birds to- 
gether with consecutive right-beak/left-beak pairs. Note that between the 
points of a component there may be points belonging to other components, 
but since two components cannot have crossing edges, any other such com- 
ponent must lie in a single interval 

Now consider the function F : M — )■ Z defined by F{0) = 0, and 

F{y)-F{x)= Yl (Lt-Rt), x<y. 

Thus, F takes a up-step (of size 2) at each left-beak and a down-step at each 
right-beak. Hence it is a continuous-time simple symmetric random walk on 
the even integers. On the other hand, by the observations above concerning 
components, if there is an infinite component, then F is bounded above a.s. 
by some (random) constant, which is impossible. □ 

Proof of Theorem \1.2\ (ii). We use a variant of an argument from [TJ Theo- 
rem 2(b)]. For A C M, write R{A) for the total number of right-stubs at 
points X G [P] n A, that is, R{A) = J2x&[v]nA^xy define L{A) analo- 
gously as the total number of left-stubs in A. Furthermore, for A, B G M., 
let R{A —J- B) denote the number of right-stubs in A that are matched with 
left-stubs in B, and D{A ^ B) the total number of edges connecting points 
in [V] n A to points in [V] fl B. Write k for the supremum of the support of 



11 



/i. For t > 0, we have 



ER{[0,2t]^[0,2tY) = l-ED{[0,2t]^[0,2tY) 



< 



2 

k '•2* 



/ F*{X > X A{2t- x)) dx 
Jo 



2 

= k-E*[XAt]. 

Furthermore, since has bounded support, we can use the central hmit 
theorem to get that 



Ei?([0, 2t] [0, 2tY) > E (i?[0, 2t] - L[0, 2t]) 



ct 



1/2 



as t — oo for some c > 0. Hence t ^^'^E*[X At] > c' for sufficiently large t and 
some c' > 0. On the other hand, if E*[X^/^] < oo, then t-^/^E*[X At] ^0 
as t — !■ oo by the dominated convergence theorem, a contradiction. □ 

Proof of Theorem M.'^ (i). Let 

H = {x^\V]: A4. > |x| - 1}, 

that is, H is the set of vertices that desire some point in the unit interval 
(— 1, 1). Write for the cardinality of H . We will show that N = oo a.s. By 
symmetry this implies that with positive probability infinitely many vertices 
in (1, oo) desire 1. However, on the latter event, for any a > 1, infinitely 
many vertices in (a, oo) desire a, so by ergodicity it follows that N = oo a.s. 

First we show that P(A^ = 0) = 0. Assume for contradiction that 
P(A^ = 0) > 0. For a configuration (V,^) with = 0, consider a modified 
configuration where a vertex is added uniformly at random in [0, 1] indepen- 
dently of V. If follows from [71 Lemma 18] and a straightforward modification 
of [TJ Lemma 16] that all stubs at this vertex would be unmatched in the 
stable multi-matching, which contradicts [Sj Proposition 2.2]. 

Now assume that all components are finite a.s., and suppose that for a 
contradiction that P(A^ < oo) > 0. For a configuration (V,^) with A^ < 
oo, consider a modified configuration where the vertices in H are removed, 
along with all vertices in their components. The number of vertices that are 
removed is almost surely finite. But in this configuration, we have A^ = 0, 
which is a contradiction to a straightforward modification of [71 Lemma 18]. 

□ 
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Proof of Theorem \1.3\ (ii). First note that, when the only possible values for 
the degrees are 1 and 2, the stable multi- matching cannot contain any cross- 
ing edges. lia<c<h<d and the edges (a, h) and (c, d) are present, then, 
as pointed out in the proof of Proposition ll.H the edge (c, h) must also be 
present in the matching. But iih — a > d — h, then h and d desire each other, 
and are hence connected by an edge, so h has degree at least 3. Similarly, if 
b — a < d — b then c has degree at least 3. 

Lemma 1231 and the fact that edges do not cross imply that an infinite com- 
ponent must consist of a set of degree-2 vertices, unbounded in both direc- 
tions, with an edge between each consecutive pair. It follows from Lemma [23] 
that the number of vertices in this infinite component that desire the origin 
is finite almost surely. 

As for the finite components, each must be contained in a single interval 
defined by an edge of the infinite component (since there are no crossing 
edges). Note also (although this observation will not be needed) that a com- 
ponent of size k must consist of vertices Xi < ■ ■ ■ < Xk with edges {xi, Xj+i) 
for alH = 1, . . . , A; — 1 together with the edge {xi,Xk)- 

Now let /q denote the interval defined by the edge in the infinite path that 
crosses the origin. This interval is finite and hence contains almost surely 
finitely many points of [P] in finite components. These points might desire 
the origin. A vertex x > (respectively, x < 0) in a finite component outside 
this interval however cannot desire the origin: if it did, it would also desire 
the left-most (right-most) end-point of the interval defined analogously to 
/q. But this vertex also desires x, which means that there would be an edge 
between them. 

We conclude that < oo almost surely, as desired. □ 

Proof of Proposition \1.5[ Say that a point x G [V] is t-bad if Mj, > t. If 
D < k almost surely, then there can be at most k t-bad points in the ball 
5(0, t/2). Hence 

k > E[number of t-bad points in 5(0, t/2)}] = vol(5(0, t/2)) P*(A/ > t), 

giving the result. □ 

Proof of Proposition Assume that the number of edges that cross the 
origin is finite with positive probability. On the event that the origin is 
crossed by finitely many edges, the same is true for any other a; G M, since 
the difference between the number of edges crossing x and the number of 
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edges crossing the origin is bounded above by the total degree of the vertices 
between the origin and x. Now consider the intervals between the points 
X G [P] with odd degrees. When passing a point with odd degree, the num- 
ber of crossing edges changes parity, that is, if it is even (odd) immediately 
to the left of the point, it is odd (even) to the right. When passing a point 
with even degree on the other hand, the parity of the number of crossing 
edges remains unchanged. This means that we can assign the value (even 
number of edges crossing) or 1 (odd number of edges crossing) to the inter- 
vals separating the odd degree vertices in a deterministic way (indeed, the 
stable multi-matching is a factor). Furthermore, the odd degree vertices con- 
stitute a Poisson process. Now, [TJ Lemma 11] asserts that it is impossible 
to assign alternating values and 1 to the intervals separating the points 
of a Poisson process as a factor of the Poisson process. Here we need the 
stronger statement that this cannot be done even using the randomness in 
the degrees of the vertices and in the position of the even degree vertices. 
This however follows from a straightforward modification of the proof of [TJ 
Lemma 11]. □ 



4 Percolation for the stable multi-matching 
with D = 2 

If the stable multi-matching almost surely has an infinite component, then 
there is a strictly positive probability p that a given vertex belongs to this 
component. Simulations of the stable multi-matching with D = 2 on large 
finite cycles indicate a largest component comprising about 0.3 of the vertices 
(see the top row of Table [Din Section |6]below). This suggests the existence of 
an infinite component with p ^ 0.3. In this section we show that percolation 
indeed follows from the assumption that a certain finite event has sufficiently 
large probability. Furthermore, we give overwhelming statistical evidence for 
this assumption. 

The key concept for the proof is the core (stable) multi-matching, 
which we define next (in the more general setting of arbitrary dimension and 
numbers of stubs). Let 5* C M*^ be a bounded set, let P C S' be a finite 
set of points, and let {Dx)x&p be positive integers representing numbers of 
stubs. Let P = P U {S^}, where := R'^ \ S. (We will treat like an 
additional point; it will not form part of the matching, but will affect the 
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notion of closest points.) Assume that all distances between pairs of elements 
of P are distinct. Assign stubs to each point x E P, and one stub to S**^. 
Repeat the following operations. From each point x G P that currently has 
an unused stub, assign an arrow pointing to the closest other element of P 
among those that have at least one unused stub and do not already have an 
edge to X. Then, for every pair x,y E P whose arrows point to each other, 
connect them with an edge and remove one stub from each. Erase all arrows 
and repeat. After some finite number of such iterations, no more edges are 
added. The core multi-matching of (P, D) in S defined to be the resulting 
graph. Note that the degree of a vertex a; G P is at most Dr^, but may be 
strictly less. 

Lemma 4.1. Let P be any discrete set of points in Mf^, let {Drc)xep positive 
integers, and let S be a bounded set. Every edge in the core multi-matching 
of (P n S, D) in S is present in every stable multi-matching of (P, D) on W^. 

Proof. It is straightforward to check by induction on the steps of the above 
algorithm that every edge added is present in any stable multi-matching. 
The key point is that if x G S* is closer to some other point in y E S than to 
S^, then X is also closer to y than to any point in P \ 5*. □ 

We now specialize to the main case of interest. Let = 1, and let = / = 
[a, 6], a bounded interval. Let P be a Poisson process of intensity 1 on M, and 
consider the case /i({2}) = 1 of deterministically two stubs per vertex. By 
the core multi-matching on the interval / we mean the core multi-matching 
of {[P] n /, -D) on /, where D = 2. We call an interval / = [a,b] good if 
the core matching on / has a connected component with a point in the first 
quarter [a, |a + ^b] and a point in the last quarter [^a + |6, b]. 

Theorem 4.2. Let V be a Poisson process of intensity 1 on M. and let 
jj,{{2}) = 1. If for some L we have P([0,L] is good) > 0.968, then the 2- 
stub stable multi-matching has an infinite component. 

Monte-Carlo simulations provide overwhelming evidence that the condi- 
tion in Theorem 14.21 indeed holds for some large L, subject to the trust- 
worthiness of the pseudo-random number generator and the software used. 
Indeed, in 1000 independent simulation runs of the process with L = 13000, 
the interval [0, L] was good in 991 cases, implying that the hypothesis that 
the probability is 0.968 or less can be rejected at the 10~^ level. See the 
appendix for details. 
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By a monotone path in a multi-matching we mean a sequence of vertices 
Xi < X2 < ■ ■ ■ < Xk with the edges (xi, X2), {x2, X3), . . . , {xk-i, Xk) aU present. 
As observed in the proof of Theorem 11.31 (ii), no two edges cross in the 
2-stub stable multi-matching, and hence the same holds in a core multi- 
matching. If / is good, it follows that the core multi-matching on / contains 
a monotone path from the first quarter to the last quarter. We call such a 
path a spanning path of the good interval. 

Lemma 4.3. Let a<b<c<dbe points ofV, and suppose that the intervals 
[a, b] and [c, d] are both longer that [b, c] . If the 2-stub stable multi-matching 
has a monotone path a from a to b and a monotone path S from c to d, then 
it has a monotone path from a to d which contains a and 5. 

Proof. Suppose on the contrary that there is no monotone path from a to d 
containing a and 5. First extend the path a to the right as far as possible 
within [6, c] ; that is, let b' G [6, c] be as large as possible such that there 
is a monotone path containing a from a to b' . Similarly extend 7 as far 
left as possible to c' G c]. By our assumption, and since there are no 
crossing edges, we have b' < d . Note also that [a, b'\ and [c', ci] are longer 
than [6',c']. Now b' is adjacent to its neighbour in the monotone path from 
a, and to exactly one other vertex x. By our assumptions, x ^ [fe',c'], and 
therefore x ^ (a, c?), otherwise we would have crossing edges. A similar 
argument shows that c' has a neighbour outside (a,(i). But now (6',c') form 
an unstable pair. □ 

Corollary 4.4. // at least 8 of the 9 intervals [0,x], [x, 2x], . . . , [8x, 9x] are 
good, then so is [0,9a;]. Furthermore, under the same assumption, given any 
spanning paths, one of each of the good short subintervals, there is a spanning 
path of the long interval containing all of them. 

Proof. Let the configuration outside I := [0, 9x] be arbitrary and consider the 
stable multi-matching. Write Ik = [{k — l)x, kx]. For any sequence of consec- 
utive good intervals /„, /q+i, ...,/&, by Lemma I4l3] we obtain a monotone path 
in their union reaching to within distance x/4 of each end. If I3, 14, 15, Iq, I7 
are all good, then the resulting path reaches to within 2x -\- x/A = 9x/4 of 
each end of /, as required. On the other hand, if one of I3, 1^, . . . , If is bad 
(but the other 8 subintervals are good), then we obtain two paths of length 
greater than 2x — 2{x/4) = 3x/2 on either side of the bad subinterval, with a 
gap of length less than x-\-2{x/4) = 3x/2 in between, so another application 
of Lemma 14.31 provides the required spanning path. □ 
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Proof of Theorem \4.2\ Let 

[ 2 ' 2 J' 

and let pk be the probability that is good. Thus po > 0.968, and by 
Corollary 14. 4^ Pk+i > fiPk) where 

/(p):=/ + 9/(l-p). 

It follows by an elementary computation that — >■ 1 as — !■ oo, and indeed 
^^(1 — pk) < oo. Hence by the Borel-Cantelli lemma, a.s. is good for all 
sufficiently large k. Moreover, since ^^(1 — fipk)) < oo, for all sufficiently 
large the interval can be divided into 9 equal intervals of which at least 
8 are good. By Corollary 14.41 it follows that for some (random) K we may 
find monotone paths t:ki'^k+1i'^k+2i ■ ■ ■■, each contained in the next, with 
TTfc a spanning path of for each k. Then Ufc>A'^fc infinite connected 
graph in the stable multi-matching. □ 



5 Some counterexamples 

Theorem 1 1 . 2 1 assert s that the random direction stable multi-matching has no 
infinite component when all vertices have degree 2, and that it has long edges 
in the sense that E*[X^/^] = oo. Furthermore, it follows from Theorem 11.31 
that existence of an infinite component in the stable-multi matching with all 
degrees equal to 2 is equivalent to < oo. This might lead one to suspect 
that there is a simple relation between the component structure and the edge 
length for /i({2}) = 1 that holds for any matching scheme. Below, we give 
two examples of factor matching schemes that demonstrate that this is not 
the case. 

Example 1. Our first example is a matching scheme where all components 
are infinite and where also the number of edges crossing the origin is a.s. 
infinite. Note that, if the origin is crossed by infinitely many edges, then 
also = oo and thus, by Lemma [1.4^ E*[M] = oo. Existence of an infinite 
component hence does not imply short edges in any of these respects. 

To describe the matching scheme, let each point in \P] be equipped with 
two stubs. Recall that the stable multi-matching in the special case where 



17 



/i({l}) = 1 is known as the stable matching. First use one stub per point 
to form edges according to the stable matching of the points. Then orient 
the remaining stub at each point in the opposite direction (left or right) 
from that of the first stub, and connect these directed stubs according to 
the procedure used for the random direction stable matching. This gives a 
graph where each point has one edge connected to the right and one edge 
connected to the left - that is, all points are birds in the in the terminology 
used in the proof of Theorem 11.21 - which implies that all components in the 
graph are infinite. That the number of edges crossing the origin is infinite 
almost surely follows from Proposition 11.61 applied to the configuration after 
the first stub per point is connected. □ 

Example 2. The next example is a matching scheme that gives almost 
surely only finite components and where the expected edge length is finite. 
Finite expected edge length hence does not imply existence of an infinite 
component. 

The matching scheme proceeds by dividing the vertices into groups of size 
at least 3 as follows. Call a point of P a seed if it has some other point within 
distance 1. Call a seed x good if the number of non-seed points between x 
and the next seed to its right is at least 2. Now whenever x < y are two 
consecutive good seeds, let all the points in [x, y) constitute one group. 

Define the matching as follows. For a group Xi < ■ ■ ■ < x^-, connect the 
two stubs per vertex to form the edges (xj,a;j+i) for i = 1, . . . , k — 1 and 
the edge This clearly gives a configuration with almost surely finite 

components and finite expected edge length. □ 

Next we give simple examples of translation-invariant point processes on 
M for which the stable multi-matching in the case /i({2}) = 1 provably does, 
and does not, have an infinite component. 

Example 3. Let (Xj)jez be i.i.d. and uniformly distributed on [0, 1/3], and 
let U be independent and uniform on [0, 1]. Consider the point process with 
support {i + Xi + U : i e Z}. It is easy to see that each point connects to 
its left-neighbour and its right-neighbour, so there is an infinite component. 

Example 4. Let {Xij)i^z,j=i,2,3 ^e i.i.d. and uniformly distributed on 
[0, 1/3], and let U be independent and uniform on [0, 1]. Consider the point 
process with support {i + Xij + U : i G Z, j = 1, 2, 3}. Then each component 
has size exactly 3. 
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Number of points 


olO 


ol2 


ol4 


9I6 




Expected degree 


2 


.244 it .099 


.291 ± .044 


.297 ± .014 


1 ^^^^^^ 

.287 ± .009 


.292 it .005 


3 


.278 ± .099 


.154 ± .029 


.049 ± .006 


.017 ± .003 


.006 ± .001 


4 


.802 ± .158 


.728 ± .232 


.653 ± .201 


.366 ± .207 


.399 ± .143 


5 


.974 ± .018 


.933 ± .114 


.815 ± .183 


.672 ± .258 


.321 ± .102 


6 


.989 ± .009 


.990 ± .004 


.975 ± .047 


.933 ± .161 


.755 ± .207 


7 


.994 ± .008 


.997 ± .003 


.989 ± .020 


.996 ± .000 


.961 ± .112 


2.1 


.071 ± .023 


.024 ± .005 


.009 ± .002 


.003 ± .000 


.001 ± .000 


2.5 


.132 ± .062 


.043 ± .018 


.018 ± .006 


.006 ± .001 


.001 ± .000 


3.5 


.472 ± .172 


.244 ± .071 


.146 ± .061 


.050 ± .014 


.016 ± .004 


4.5 


.992 ± .011 


.888 ± .110 


.529 ± .138 


.298 ± .132 


.129 ± .059 



Table 1: Simulation results for the stable multi-matching of uniformly ran- 
dom points on the cycle. The proportion of points in the largest connected 
component is indicated as "sample mean ± sample standard deviation" for 
a sample of size 10. The degree D is either a constant integer, or takes two 
consecutive integer values with probabilities determined by the indicated ex- 
pected value. 

6 Open problems 

The random direction stable mult i- matching 

For degree distributions other than yU,({2}) = 1, it remains an open problem 
to determine if the random direction stable multi-matching generates an 
infinite component. 

The stable mult i- matching 

Firstly, it would of course be desirable to turn the "statistical proof" of perco- 
lation for D = 1'm.d = \ into a fully rigorous proof. Furthermore, it remains 
an open problem to determine whether there exists an infinite component 
the stable multi-matching for other degree distributions (an exception being 
the case D G {1,2} with P(D = 1) > 0; see [5]). Another interesting case 
arises when most vertices have 2 stubs, but a small fraction have 3; this case 
can be expected to be very different from the 2-stub case since there are lo- 
cal configurations which can end a long path. Indeed, simulations appear to 
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indicate that the proportion of vertices in the largest component converges 
to as the system size increases, thus suggesting no infinite component; see 
Table [1] (fines 9-10). Tfie observation tfiat local modifications can end a long 
path extends to any degree distribution with support on at least one odd 
integer. Is it the case that the stable multi-matching has an infinite compo- 
nent in d = 1 if and only if the degree distribution has support only on even 
integers? The results in Table [T] appear consistent with this hypothesis. 

Iterated stable matching 

Yet another mult i- matching scheme is obtained by repeatedly applying the 
stable matching of the points with the restriction that multiple edges are 
not allowed. More specifically, first take the stable matching of [P], connect 
the points accordingly and remove one stub per point. Then consider the 
stable matching of the points that have at least one stub left on them and 
with the modification that two points that already have an edge between 
them cannot be matched. This matching is obtained by repeatedly matching 
mutually nearest neighbors in the set of points with at least one stub left on 
them, avoiding matchings of points that are already connected. As remarked 
in [21 Remark 2.2], the proof of [21 Proposition 2.2] is easily modified to show 
that this yields a perfect matching. Connect the points according to this 
matching and remove one stub from each point that is connected. Repeat 
indefinitely. 

Does this matching scheme generate an infinite component? How does 
the answer depend on the degree distribution? Note that it follows from 
Proposition 11.61 that the number of edges that cross the origin is a.s. infinite 
already after the first stub of the vertices has been connected. For degree 
distributions with degrees larger than 1 however the matching contains cross- 
ing edges. This means for instance that the proof of Theorem 11.3( b) cannot 
be applied to draw the same conclusion (that < oo if there is an infinite 
component) for the iterated stable matching. 
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Appendix — simulation code 

The Python 2.6 code below was used to verify that the interval [0, 13000] was 
a good interval in 991 out of 1000 cases. Since 

P[Binomial(1000, 0.968) > 99l] < 10"^ 

this gives grounds to reject the hypothesis that P([0, 13000] is good) < 0.968 
at the 10-6 level. 

The code uses Python's built-in implementation of the Mersenne Twister, 
one of the most extensively tested pseudo-random number generators. The 
experiment was also repeated using two other pseudo-random number gener- 
ators, and using an alternative method of generating a Poisson point process, 
with consistent results. 

The following observations simplify and speed up the construction of the 
core multi-matching in the 2-stub 1-dimensional case. In place of the com- 
plement of an interval [a, b], it suffices to consider distances to the endpoints 
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{a, b}. At any stage of the core multi-matching algorithm, the arrow from a 
point X with an unused stub must point to one of x's two closest neighbours 
on the left or its two closest neighbours on the right among the set of points 
of P with unused stubs; this is because at most one other point can have an 
edge to X already. Also, if the arrows of x, y point to each other, then there 
is no other point z with unused stubs between them; indeed such a z cannot 
already have edges to both x and y, so one of x and y would instead point 
to z. 



from random import * 
from math import * 



def poi(a): # Poisson random variable with mean a 

t=-l 

while a>0 : 

a+=log(random() ) 

t+=l 
return t 



def setup (a) : 

global n,x, stubs, e 
n=poi (a) 

x= [random () for i in xrange(n)] 

x.sortO # X = sorted list of n=poi(a) random points in [0,1] 

stubs=[2]*n # 2 stubs per point 
e=set([]) # e = set of edges 

def xx(i) : # position of point i, or endpoint of [0,1] for i outside range 

if i<0: 

return . 
elif i>=n: 

return 1 . 
else : 

return x[i] 

def corematchO : 
cont=True 
while cont : 

cont=False 

active= [-1 ,-1] + [i for i in xraiige(n) if stubs [i] ] + [n,n] 

# points with unused stubs, plus 2 dummy points at each end 
arrow= [None] *len(active) 
for j in xrange (2 , len(active) -2) : 

1= [Cabs (xx(active [j] ) -xxCactive [k] ) ) ,k) # find distances 

for k in j-2 , j-1 , j+1 , j+2 # to 2 neighbours on each side 

if tuple (sorted( [active [j] , active [k] ]) ) not in e] # if no edge already 
if 1: 

arrow[j]=min(l) [1] # arrow to closest 

for j in xrange (2 , len(active) -3) : 

if arrow [j] ==j+l and arrow [j +1] ==j : # found pair with mutual arrows 

# (must be neighbours) 
e . add( (active [j] , active [j+1] ) ) # add edge and remove stubs 

stubs [active [j] ] -=1 
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stubs [active [j+1] ] -=1 
cont=True 



# keep going if some edge added 



def components : 

nbrs=dict ( (i , [] ) for i in xrange(n)) 
for (i,j) in e : 

nbrs [i] .append(j) 

nbrsEj] . append (i) 
done= [False] *n 
ans= [] 

for i in xrange(n): 
if not done [i] : 
cur= [i] 
done [i] =True 
for j in cur: 

for k in nbrs[j] : 
if not done [k] : 
cur . append (k) 
done [k] =True 
ans . append(sorted(cur) ) 
return ans 



# find components of graph 



def goodO : # is the interval good? 

return any (x [c [0] ] < . 25 and x[c[-l]]>.75 for c in componentsO) 

seed(12345) 
a=13000 
k=1000 
g=0 

for i in xrange(k): 
setup(a) 
corematchO 
if goodO : 
g+=l 

print g, '/' , (i+1) 
print 

print 'Interval of length' , a, 'was good' ,g, 'times out of',k 
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